Hyperbolicity of scalar-tensor theories of gravity 
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Two first order strongly hyperbolic formulations of scalar-tensor theories of gravity allowing non- 
minimal couplings (Jordan frame) are presented along the lines of the 3+1 decomposition of space- 
time. One is based on the Bona-Masso formulation, while the other one employs a conformal 
decomposition similar to that of Baumgarte-Shapiro-Shibata-Nakamura. A modified Bona-Masso 
slicing condition adapted to the scalar-tensor theory is proposed for the analysis. This study con- 
firms that the scalar-tensor theory has a well posed Cauchy problem even when formulated in the 
Jordan frame. 
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I. INTRODUCTION 

Scalar-tensor theories of gravity (STT) are alternative 
theories of gravitation where a scalar field is coupled non- 
minimally to the curvature associated with the physical 
metric (this is the so-called Jordan frame representa- 
tion). The term "physical metric" refers to a situation 
where test particles follow the geodesies of that metric. 
The variation of the action of the STT with respect to 
the physical metric gives rise to field equations which 
contain an effective energy-momentum tensor (EMT) in- 
volving second order derivatives in time and space of the 
scalar field. Such EMT has the property that "ordinary 
matter", i.e. matter associated with fields other that 
the scalar field, obeys the (weak) equivalence principle 
which mathematically translates into a conserved EMT 
for ordinary matter alone. 

Since a priori it was not clear how such second or- 
der derivatives could be eliminated (in terms of lower 
order derivatives) , or managed so as to obtain a quasilin- 
ear system of hyperbolic equations for which the Cauchy 
problem was well-posed (in the Hadamard sense), many 
people decided to abandon this approach in favor of the 
so-called Einstein frame representation where the non- 
minimal coupling (NMC) is absorbed into the curvature 
by means of a conformal transformation of the metric. 
The new conformal metric is unphysical in the sense 
that (non-null) test particles do not follow the geodesies 
of that metric. However, the mathematical advantage 
is that the field equations for the nonphysical metric 
resemble the standard Einstein field equations with an 
unphysical effective EMT which involves at most first 
order derivatives of a suitable transformed scalar field 
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(this EMT is unphysical because the "ordinary matter" 
part is not separately conserved). In the Einstein frame 
one can show that by using standard gauges {e.g. har- 
monic gauge) the field equations acquire the form re- 
quired in the application of theorems that establish the 
well-posedness of the Cauchy problem. 

In view of the apparent mathematical advantages and 
disadvantages of the Jordan and Einstein frames, sev- 
eral widespread misconceptions became common in the 
literature. One of these concerned the statement that 
the Cauchy problem is only well-posed in the Einstein 
frame [l|. In Ref. Q, however, one of us showed that 
this is not the case by following two different approaches. 
One was in the spirit of a second order analysis consisting 
on reducing the set of field equations, both for the metric 
components and the scalar field, to a manifestly quasilin- 
ear diagonal second order hyperbolic form (the "reduced" 
field equations). This was achieved by manipulating the 
field equations in a way that allowed one to express the 
d'Alambertian on the scalar field in terms of at most first 
order derivatives, and also by implementing a modified 
harmonic gauge which was adapted to the STT. Such 
a modified gauge allowed one to eliminate the remain- 
ing second order derivatives of the scalar field which had 
previously prevented the applicability of Leray's theorem 
(see e.g. Ref. @ for the theorem). 

The second approach followed in Q, and which was 
directly related with the initial value problem, consisted 
in recasting the field equations of the STT in a 3+1 
or Arnowitt-Deser-Misner (ADM) [4| form written a la 
York Q (hereafter referred to as the ADMY equations). 
To do so it was necessary to define suitable variables, fol- 
lowed by a manipulation of the equations in order to ob- 
tain well-defined constraints (independent of the choice of 
the lapse and shift), as well as first order in time evolution 
equations for the extrinsic curvature and the scalar-field 
variables. In Ref. Q the main goal was to focus in the 
3+1 approach rather than the second order one because 
of the subsequent numerical applications we had in mind. 
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In fact, almost all the modern codes used in numerical 
relativity are based on first order in time formulations, 
and therefore we wanted to adapt them for the analysis 
of phenomena in STT. 

Now, when taking the limit of pure general relativity 
(GR) (i.e. no no-minimal coupling), the 3+1 equations 
presented in reduce to the standard ADMY equations 
with a minimally coupled scalar field plus ordinary mat- 
ter sources. It is well-known that the ADMY equations 
are not strongly hyperbolic (see [f| for a detailed discus- 
sion), and therefore the corresponding equations for STT 
described in Q were expected not to be strongly hyper- 
bolic either. As already mentioned in Q, such equations 
were only the first step towards a first order strongly 
hyperbolic system for which the well-posedness of the 
Cauchy problem could be established. 

The aim of this paper is then to fill that gap, and to 
obtain a first order strongly hyperbolic system of par- 
tial differential equations based on the 3+1 equations 
of 0]. Actually, we will show here two such systems: one 
that is related to the Bona-Masso-Seidel-Stela (BMSS) 
formulation 0] (which in turn is based on the Bona- 
Masso (BM) formulation |8|), and a second one based on 
the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) con- 
formal decomposition 0, [l(| ■ As was done in [|J , a modi- 
fied BM slicing condition will be used for the STT, while 
the shift vector will be taken as an a priori given function 
of the coordinates. 

A last comment about the Jordan vs. Einstein frame is 
in order. In [2j it was a matter of principle to show that 
the Cauchy problem could be well formulated in the Jor- 
dan frame. It then became clear that there was no fun- 
damental reason to continue using the Einstein frame. 
After all, the Jordan frame is the one associated with 
the (physical) quantities that are to be confronted with 
the observations. Moreover, the constraint and evolution 
equations are not more involved than those of usual GR, 
so from a numerical point of view the use of the Jordan 
frame does not add much complexity to the analysis. On 
the other hand, the direct use of the Jordan frame al- 
lows a better physical interpretation of the results and 
avoids the potential problems that might arise in the 
Einstein frame in cases when the conformal transforma- 
tions back to the physical metric are not well-defined. 
Finally, it has been argued that some of the energy den- 
sities associated with the Jordan frame effective EMT 
(JFEMT) are not positive definite and may therefore be 
unphysical, and also the corresponding ADM mass could 
then become negative. While it is certainly true that the 
JFEMT does not satisfy the energy conditions in general, 
it is very likely that any physical configuration consistent 
with observations will not carry (total) negative energy. 
This is because the deviations of any alternative theory 
from GR are presumably very small in order to repro- 
duce many of the current observations (cf. Ref. [111]). 



Actually, in any viable STT cosmology it turns out that 
the total energy density of the Universe is always posi- 
tive, even when taking into account the negative contri- 
butions due to the NMC [H, [H, [H, \tf& . Again, this is 
because the scalar-field contributions (positive or nega- 
tive) to the total energy density should be in agreement 
with several observations related to the past and present 
history of the Universe. For instance, the expansion rate 
of the Universe modified by the NMC contributions has 
to be consistent with the one required to produce the 
correct abundance of primordial nucleosynthesis [l2T |. In 
addition, the density perturbations in STT should also 
match the Cosmic Microwave Background data. In other 
words, a STT cosmology which produces a total effec- 
tive negative energy density of the Universe at any given 
epoch will surely not be consistent with observations. By 
the same token, in a consistent STT cosmology the nega- 
tive contributions (if any) are to be naturally suppressed 
by the positive ones, providing a consistent Universe. 

In the case of astrophysical applications (e.g. neutron 
star models) within a class of STT with a positive defi- 
nite NMC function (which implies a positive definite ef- 
fective gravitational "constant"), the ADM mass turns to 
be always positive despite the negative contributions to 
the energy density due to the NMC [Hj]. Moreover, the 
value of the NMC cannot be very high as otherwise the 
corresponding STT would put in jeopardy the agreement 
with the binary pulsar observations [l7| . 

Of course, the nice thing about matter fields satisfying 
the energy conditions is the applicability of several the- 
orems (e.g. positive mass theorems and singularity the- 
orems; see Q and references therein). However, the fact 
that the JFEMT does not satisfy the energy conditions 
in general simply indicates that the theorems cannot say 
anything about the positivity of mass or the formation of 
singularities in this case. This is not a problem of physical 
but of mathematical character. But again, the nonposi- 
tive definite terms arising from the NMC will be surely 
bounded if the STT in hand is to be consistent with the 
current observations and experiments, and therefore the 
nice features of an EMT respecting the energy condi- 
tions will also very likely appear in any observationally 
consistent STT. Of course, all these arguments are only 
sustained by numerical experiments when constructing 
viable phcnomcnological models (in cosmology and com- 
pact objects), and therefore do not constitute a theorem. 
It would then be quite interesting to explore the possibil- 
ity of proving positive energy and singularity theorems in 
the Jordan frame if one can bound the non-positive def- 
inite contribution associated with the JFEMT. 

This paper is organized as follows. In Sec. [TT] we in- 
troduce the STT and the 3+1 equations described in 0]. 
The notation of several variables will be slightly modified 
relative to Ref. [2j in order to match with the one usually 
employed in numerical relativity. In Sec. IIIII a fully first 
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order strongly hyperbolic system of the field equations of 
STT is obtained along the lines of the BMSS and BSSN 
formulations of GR. We conclude with a discussion of 
future numerical applications of the hyperbolic systems 
presented here. Finally, in Appendix A we include the 
complete set of equations for both formulations, and in 
Appendix B we introduce a simple example to comple- 
ment the ideas of Sec. 



II. SCALAR-TENSOR THEORIES OF GRAVITY 

The general action for STT with a single scalar field is 
given by 



S 
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(2.1) 



with <f> the nonminimally coupled scalar field, and where 
ij) represents collectively the matter fields, i.e. fields 
other than (f>, and Go is the usual gravitational constant, 
(we use units such that c = 1). 

The representation of the STT given by Eq. ()2.ip is 
called the Jordan frame representation. The field equa- 
tions obtained from the action (|2.ip are given by [37| 
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where ' indicates d^, □ := £/ ,afc V a Vb is the covariant 
dAlambertian operator, G a b = Rab — \dabR-, and 
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Using Eq. (|2.2p . the Ricci scalar can be expressed in 
terms of the energy- momentum tensor Eq. (|2.4|) . Equa- 
tion (|2.3p then takes the following form, 



□ c 



fV - 2f'V - \f (1 + 3/") (V0) 2 + i/'Tmatt 



(2.8) 

where T m att stands for the trace of T^ b att (the subscript 
"matt" refers to matter fields other that <fi). 
Now, the Bianchi identities imply that 



V r T ca = 



(2.9) 



However, the use of the field equations leads, as men- 
tioned before, to the conservation of the EMT of the 
matter alone 



V T € 



(2.10) 



which implies the fulfillment of the (weak) equivalence 
principle. 

In what follows we shall use the 3+1 formalism of 
spacetime [1, [H[ in order to recast the field equations 
as a Cauchy initial value problem. The following quan- 
tities turn out to be useful in the 3+1 decomposition of 
the field equations of STT 0: 



3a := D a 4 
n := £n< 



(2.11) 
(2.12) 



where Cn stands for the Lie derivative along the nor- 
mal n a to the spacelike hypersurfaces St , and the opera- 
tor D a is the covariant derivative compatible with the 
3-metric "f a b induced on E t . In components we have 
n a = (l/a,—(3 l /a), where a and (5 % are the lapse func- 
tion and shift vector, respectively, associated with the 
spacetime coordinates (t,x l ) :38|. 

From definitions (|2.1ip and (|2.12p it is possible to write 
an evolution equation for Q a : 



C n Qa = -D a (aU) 
a 



(2.13) 



On the other hand, the definition of the extrinsic curva- 
ture is given by 



(2.14) 



K a b ■— — - C-IXlab 



In this way, the 3+1 decomposition of equations 
leads to a set of constraint equations plus the evolution 
equations for the extrinsic curvature (see Ref. [2| for de- 
tails). The Hamiltonian and momentum constraints are 
respectively 
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where we have denoted Q 2 — QiQ 1 , E n 
and ./matt : = ~ n dl a c T nLu and where 7 a c = 6 a c + n a n c 
is the projector onto S t . The evolution equation for the 
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extrinsic curvature takes the form: 



d t K l , - f3 l diK l } - K\d 3 (5 l + K l foP + D l D 3 a 
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where S^ b att := l a c l b d T^ att and S matt is its trace. 

Finally, equation (|2.8[) can be written as the following 
first order evolution equation (2j: 

£ n n -nif- Q c L> c [lna] - L> C Q C 

fV - 2f'V - if (1 + 3/") (q 2 - n 2 ) + |/T matt 



(2.18) 



where T mat t = 5 ma tt - -Ematt ■ Here one must remember 
that indices of 3-tensors (including 3-vectors) are raised 
and lowered with the 3- metric 7, and that the contravari- 
ant time components of such quantities are identically 
null. This is why one is usually only interested in the 
spatial components of 3-tensors. 

A useful evolution equation for the trace of the extrin- 
sic curvature takes the form: 



d t K - (3 l diK + 3 Aa - aK io K i3 - (DiQ 1 + UK 




(i+2r)-V 
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(2.19) 



where 3 A := D l Di is the Laplacian compatible with the 
3-metric, and D C Q C = d c Q c + Q c d c (ln^f) (here 7 = 
det/yy). 



The system of equations ([2TT2"]) . (|2TT5|l and (|2~T4"]) can 

now be rewritten as 

9t<t> = aU + (3 a Q a , (2.20) 
d t Q l - (3 l diQi - Qidtf 1 = A(oH) , (2.21) 
dtHj = -2(aK l3 - D {i p j} ) . (2.22) 

This system of evolution equations is to be completed 
with appropriate evolution equation for the gauge vari- 
ables (lapse and shift). This issue is considered next. 

In order to obtain a closed evolution system one has 
to impose gauge conditions for the time variable t and 
for the spatial coordinates x l . In this work we shall con- 
sider a modified Bona-Masso time slicing defined in the 
following way (cf. Ref. (H) 



Ut = 



1 



BM 



- 1 n"n c V h V r t 



e f 

JBM J 



(2.23) 

with Jbm — Ibm{o) > a positive but otherwise arbi- 
trary function of the lapse, and O = 8(/bm) is in princi- 
ple an arbitrary function of Jbm which can a posteriori 
be fixed so as to ensure a well behaved hyperbolic system 
(see Section|TlT]below). It is important not to confuse the 
gauge function /bm(o) with the NMC f(4>). 

Concerning the spatial coordinates, we shall consider 
the shift vector as an a priori known function of the coor- 
dinates. However, in the future it would be interesting to 
investigate some "live" shift conditions and their effects 
in phenomena involving STT. 

Now, the slicing condition (|2.23p for local coordinates 
x a = (t, x l ) adapted to the 3+1 foliation of the spacetime 
reduces to 



da 2 



BM 



K 



IBM f 



(2.24) 



where d/dt := dt — £p. The specific choices = Jbm 
1 correspond to a modified harmonic slicing condition 
(termed "pseudoharmonic" in (2j), which was specially 
useful for the second order hyperbolicity analysis per- 
formed in (2j. On the other hand, with = one recov- 
ers the usual BM slicing condition. Just like in standard 
GR, we shall see that one can take in general any posi- 
tive Jbm (not necessarily Jbm = 1) to obtain a well be- 
haved strongly hyperbolic system. Moreover, the choice 
= 1 or = fsM provides simple sufficient conditions 
for regular eigenfields in the STT, as we will show be- 
low. An inadequate choice for leads to eigenfields that 
can be ill-defined when Jbm = 1 locally or globally (cf. 
Eqs. (|3~24|) and (f3T49|) below). A bad choice can be pre- 
cisely = 0, so the usual BM slicing condition seems not 
to be well adapted for constructing strongly hyperbolic 
formulations of STT's in the Jordan frame. 
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III. HYPERBOLIC SYSTEMS 

In the previous Section we showed a constrained sys- 
tem of first order in time partial differential equations 
(PDE) corresponding to the STT. This system is well 
formulated in the sense that one can construct an unam- 
biguous numerical algorithm that evolves initial data sat- 
isfying the constraint equations. In the case of a globally 
null scalar field, the above system of equations reduces 
to the usual ADMY system of equations of GR. Even if 
such a system is well formulated in the above sense, it 
is well-known that the evolution is not stable in many 
cases. 

It took several years to realize that such instabilities 
were not associated with the numerical algorithms but 
rather with the mathematical structure of the ADMY 
system itself. During the mid 90's, many researchers 
started to suspect that the numerical instabilities that 
seemed to plague several testbed codes based in the usual 
3+1 formalism were caused by the lack of (strong) hyper- 
bolicity of the system. Therefore, many investigations 
were devoted to develop new strongly hyperbolic formu- 
lations of the 3+1 evolution equations of GR (see [9| and 
references therein). However, at that time, it was still 
not totally clear that such new hyperbolic formulations 
presented more computational advantages than the usual 
ADMY system. The breakthrough towards a real im- 
provement of the numerical stability was the use of a 
novel formulation of GR which used a conformal formal- 
ism, today known as the BSSN formulation [1, [l(|. In Q 
it was recognized that the new conformal system pre- 
sented some hyperbolic features not present in the orig- 
inal ADMY system, and moreover it was stressed that 
the use of suitably defined connection coefficients and 
the subsequent addition of the momentum constraints 
(in order to eliminate some divergence operators) proved 
to be essential for numerical stability. Nonetheless, it 
was after a more systematic analysis of the BSSN system 
that a much better understanding of its hyperbolicity 
properties was obtained [19( ] . In connection with the hy- 
perbolicity properties of the 3+1 system, one needs to 
mention the gauge choice (i.e., the prescription of the 
evolution of the gauge variables) . In the formal develop- 
ments and proofs of theorems about the well-posedness 
of the Cauchy problem of GR, it was the use of the har- 
monic gauge that allowed one to write the Einstein field 
equations as a second order diagonal hyperbolic system 
of quasilinear PDE's. Nevertheless, in numerical imple- 
mentations, the use of the (pure) harmonic gauge has 
proven not to be very successful, specially in the strong 
gravity regime where one requires to avoid (or freeze) 
the formation of singularities in the numerical grid. A 
big step towards solving this problem was the imple- 
mentation of new gauge conditions (e.g. the so-called 
Bona-Masso time slicing) in hyperbolic formulations 



Currently, many hyperbolic formulations, together with 
Bona-Masso like slicing conditions, have been proposed 
in the literature. 

In this and the following sections we will propose two 
systems for the STT based on the the 3+1 system of 
Sec. |TT] but written in fully first order form along the 
lines of the BSSN and the BMSS formulations, and show 
that they are indeed strongly hyperbolic. Our aim for 
constructing two strongly hyperbolic formulations of the 
STT is to understand the way in which the modified 
Bona-Masso condition (|2.24[) enters into both formula- 
tions, and how robust is this new class of slicing condi- 
tions. In fact we will see that both formulations are com- 
patible with the choices = 1,/bm m that the eigen- 
fields turn out to be well behaved without imposing any 
further stringent conditions on the gauge function fsM 
(apart from it being positive definite). 

Let us first outline the general approach we intend to 
follow for both formulations. From the numerical point 
of view, it seems to be more useful to recast a second 
order system of PDE of hyperbolic type into a full first 
order form as follows 

d t u + M%u = S(u) , (3.1) 

where u represents collectively the fundamental variables 
(like the 7y 's, -Ky's, etc.), MP are called the character- 
istic matrices of the system, and S(u) represents source 
terms which include only the fundamental variables but 
not their derivatives. This system of PDE is said to be 
quasilinear since it is linear in the derivatives but in gen- 
eral is non-linear in the u. 

A quasilinear system of PDE of the form (|3.1[) possesses 
a well-posed Cauchy problem (in the Hadamard sense) 
only in very special circumstances (see Refs. [ISHIl for a 
review), which depend on the structure of the character- 
istic matrices M* . In general, when one deals with such a 
system one can choose a specific spatial direction of prop- 
agation of fields along which one analyzes the hyperbolic 
character of the system. For instance, if the direction 
of propagation is specified by a unit covector Vi, then 
the hyperbolic nature of the system will be determined 
by the matrix C := t^M 1 , known as the principal symbol. 
Namely, the system will be weakly hyperbolic or strongly 
hyperbolic, if the following properties are satisfied respec- 
tively: 1) C has real eigenvalues but an incomplete set 
of eigenvectors; 2) C has real eigenvalues and a complete 
set of eigenvectors (for any Wj). In addition, the system is 
said to be symmetric or symmetrizable hyperbolic, res- 
pectively, if: 3) MP are symmetric, 4) MP can be sym- 
metrized. Only strongly, symmetric and symmetrizable 
hyperbolic systems of PDE admit a well-posed Cauchy 
problem. 

The systematic approach to analyze the hyperbolicity 
of a system of evolution equations along these lines con- 
sists then in first recasting the system in the form (|3.1[) . 
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and then solving the eigenvalue problem of the corre- 
sponding principal symbol. If C satisfies the conditions 
'2)' above then one can ensure well-posedness. In gen- 
eral, one finds that for the systems of evolution equations 
found in the 3+1 formulation of GR, it is necessary to add 
multiples of the constraint equations (|2. 15[) and (|2.16p 
to some of the evolution equations in order to make the 
system strongly hyperbolic [22I ] . This procedure does not 
affect the physical solutions (i.e. those that satisfy the 
constraint equations), but it does affect the structure of 
the C matrix. 

Now, suppose that the following system is strongly hy- 
perbolic in a given direction Vi (i.e., in what follows we 
will assume that C satisfies the condition '2)' above). 
From Eq. (|3.1[) one finds that 



dt.u + CS7 v il ~ 



(3.2) 



where V„ stands for the directional derivative along Vi, 
and ~ indicates that we shall focus only in the principal 
part of (|3.1j) . If we define now R as the matrix of eigen- 
vectors of C, then L = R _1 CR is the diagonal eigenvalue 
matrix. Notice that R can clearly be inverted since we 
have assumed that C has a complete set of eigenvectors. 
Multiplication of (f3~2")l on the left by R" 1 leads to the 
equivalent system 



d t w + LV s w ~ 



(3.3) 



where we have introduced the functions w := R _1 u which 
are called the eigenfields or eigenf unctions. Clearly the 
new system (I3.3|) is now decoupled and each eigenfield Wi 
propagates with its own characteristic speed given by its 
corresponding eigenvalue A^. 

All this construction depends heavily on the existence 
of a complete set of eigenvectors (in addition to the exis- 
tence of real eigenvalues). The system (|3.2p is therefore 
equivalent to (|3.3|) when the strongly hyperbolicity con- 
ditions are met. Furthermore, conditions on the smooth- 
ness of the eigenvectors (or eigenfields) and several other 
conditions concerning the gauge speeds (e.g., real- valued 
speeds, singularity avoidance) impose in turn conditions 
on the class of gauges adopted in the selection of the time 
slicing (see below) and the spatial coordinates (although 
here we shall consider the shift condition as a prescribed 
one). 

Now, it turns out that very often (depending on the 
complexity of the system) one can construct the equiv- 
alent system (|3.3p not by using direct methods to dia- 
gonalize the matrices MP, but instead by a "judicious 
guessing" approach which consists in constructing the 
eigenfields and corresponding eigenvalues by inspection 
(see Appendix B for a simple example). This is for ins- 
tance the approach followed in Q ■ In sections IIIA and 
IIIB below we use this inspection approach to construct 
the eigenfields and their corresponding speeds of propa- 
gation. 



A. Bona-Masso Formulation of STT 

In order to construct a full first order system based on 
the system of equations of Sec. II, we first need to define 
new first order variables in the following way 



a, = dj lna 



dkij = - dk lij 



(3.4) 



Moreover, we will also introduce the auxiliary variable 

Vi = di m d m i , (3-5) 

where d^ := j lm dij m and d l ij :— "f lm d m ij. Using this 
definition, the principal part of 3 Rij reads 



Ri 



— d m d m ij — d(j (ZVj) — dj) m 



(3.6) 



With such definitions the evolution equations (|2.17|) will 
be of first order in space for and dkij , instead of second 
order for the lapse and the 3-metric. 

After some algebra one obtains the following set of 34 
first order evolution equations, written here only up to 
principal part (the full form of the equations can be seen 
in the Appendix A): 



d Qi 
doll 

d CLi 

dodijk 
d Kij 

d V i 



a d l (5\ n) , 

a diQ 1 , 
-ad, xl 



-a diKjk , 
-a diA l ij , 



a di 



(1 + (K il _/if)+ c /^n 



(3.7) 
(3.8) 

(3.9) 

(3.10) 
(3.11) 

(3.12) 



where do := d t — (3 l di, V 1 :— j ll Vi and 



A z — d l ■ 



f 

jQi) 



2V, 



3) 



,( 3 - 13 ) 

Though the A y are not independent variables, it is nev- 
ertheless useful to have their evolution equation which 
take the form 



d A 1 



- ai lm d m K l3 + aS l {i d n 



2(1 + ( K 



0) 



-8 m 3) K) + 6") } (1 - f BM )K + L m m 3) (l + 6 + 20 

(3.14) 



Here we must stress the fact that in Eq. (|3 . 1 2[) we have 
added a multiple as of the momentum constraint (|2.16[) . 
Moreover, the evolution equations for the lapse and the 
3-metric were not included since they are "trivial" in the 
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sense that the lapse and the 3-metric do not evolve up to 
principal part and only they propagate along the "normal 
lines". As pointed out in [23(], the above evolution sys- 
tem, Eqs. (|3.7p — (I3.12[) . can be seen as a reduced system 
evolving in a inhomogeneous "background" . 

Now, in order to obtain the eigenfields propagating 
in one specific direction we will consider the l x' direc- 
tion, i.e. Vi — (1,0,0), and ignore the derivatives in the 
other directions (cf. Q). This is equivalent to analyzing 
only the characteristic matrix MP . Thus, for q ^ x, it 
is clear from the evolution equations (|3.7|) — (|3.12[) that 
there are 18 eigenfunctions which propagate with speed 
—f3 x , namely 101,2 = Qq (corresponding to the two eigen- 
fields propagating in the two directions orthogonal to 'x'), 
W 3A = a q, ^5_i6 = d q jk (12 eigenfields: six for each q), 
and iui7,i8 = V q + (1 + s)d l lq . 

Furthermore, by using the evolution equations we ap- 
preciate that the following combinations 

/' 

w w = a x - f B M d xm - — Q x , (3.15) 

w 20 = V x -<; f jQ x -(l + <;)A x , (3.16) 

provide two more eigenfunctions propagating along the 
normal lines with — (3 X speed, where 



A x := A x 



^XXX 



r 



--( 



A x 
p' 



, (3.17) 



on the surface x = const, (for 
Here we used A x := ^ lm A x i m 
Such a projection is in fact 



is the projection of A x pq 
q 7^ x and p / 1) [6j. 
and A xxx := -f xl j xm A x i m . 
achieved with (P x )ij := 7y — (s x )i(s x )j, where s x is the 
normal to the surface with components given by (s x )i — 
Sf/Vl^ and {s x Y = 7 r r 7 % /r er . So A x = [I")' 1 A',,. In 
fact by projecting Eq. (|3.14p on the surface x = const., 
and focusing in propagation only in the 'x' direction, 
one obtains doA x ~ —aj xx d x K , where K := (P x ) l: > Kij 
whose expression is given below in Eq. (|3.25p . 

Raising the index in Eq. (|3.7[) . and combining the evo- 
lution equation for Q x with Eq. (|3.8[) . we obtain the fol- 
lowing two eigenfunctions 



1021,22 



(3.18) 



which propagate along the light cones with eigenvalues 
A± = ~P X T a^f* ■ (3.19) 

Next, combining the equations for K qp and A x qp , we form 
the following six more eigenfunctions: 



1023-28 = A X pq ± Vl^K, 



L P9 ' 



(3.20) 



which propagate also along the light cones with eigenva- 
lues 



<j± = -f3 x ± a v 7 !^ 



(3.21) 



Combining the expressions for K x q and A xx q (A xx q 
r ) xl A x i q ), one can see that 



are four more eigenfunctions with eigenvalues 

lj± = -f3 x ± a y/-<;-f xx . 
Finally the last 2 eigenfunctions turn out to be 



_ u . 2(1 + » 

^33,34 — A + -A 

IBM — 1 

2(1 + 



(3.22) 
(3.23) 



±Vf. 



BM 7 



K 



f. 



BM 



1 



K 



2^ + i + e\ /' 



Jbm — 1 



J [Q X + VfBMT 



n 



(3.24) 



where 



K := K — 



K xa 



iB. , (3-25) 



is defined in a similar fashion as A x (cf. Eq. (|3.17p V 

The eigenfields 1033,34 propagate with the "gauge 
speeds" given by the following eigenvalues 



(3.26) 



We have then 34 linearly independent eigenfunctions 
Wi-34 which are equivalent to the 34 linearly indepen- 
dent variables (Qi, II, <2j, dijk, K^, V 1 ), so the system is 
strongly hyperbolic. 

From the 34 eigenfunctions obtained above, we appre- 
ciate that only 1^33,34 might develop gauge divergences, 
while the remaining eigenfields are smooth (as long as 
the original variables themselves are smooth) regardless 
of the time slicing. Note, for instance, how in Eq. (|3.24|) 
the choice /bm = 1 can cause those eigenfunctions to 
blowup unless we choose the value <; = — 1 (see the co- 
efficients of A x and K in 1033,34). Moreover, even with 
<T = — 1 , another blowup of the same sort might appear in 
the coefficient with y of 1033,34. However, that problem 
can in turn be avoided by a suitable choice for 8. The 
simplest choices are = 1 and = Jbm- For instance, 
the eigenfunctions 1033,34 with <; = — 1 and = 1 simply 
become 



W33 ,34 = A* ± y/f BM ~i xx K . (3.27) 
Moreover, with ? = — 1 the eigenficld (13. 16[) becomes 

/' 



1020 



f 



(3.28) 



The value s = — 1 also turns out to be compatible 
with the condition for real-valued eigenfunctions W29-32 
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(cf. Eq. (|3.22p ). Now, when one chooses fsu = 1, 
the two optimal values for O reduce both to = 1, 
and the slicing condition (|2.24[) is compatible with the 
modified harmonic time slicing which together with the 
modified harmonic spatial coordinates (both written as 
□x a = — j V/) were the ones required in Ref. [2j to 
reduce the field equations of STT (without matter) into 
a quasilinear diagonal second order hyperbolic system. 
Furthermore, note that 1^33,34 degenerate (1033 = 11134) if 
/bm = 0, and the eigenfields are not complete in that 
case. We then appreciate the importance of the condi- 
tion fsm > 0. In particular, in the BMSS formulation of 
standard GR, the condition Jbm > 1 seemed to be the 
most adequate for a convenient mode propagation and 
singularity avoidance behavior. 

Finally, let us remark that in pure GR the eigenfield 
(|3.28|) is given by w^q = V x which propagates along the 
normal lines. We just saw that u>i7,i8 = V q (with q = 
— 1) also propagate along the normal lines. Therefore 
in pure GR all three Vi's propagate along the normal 
lines. In pure GR there are also only 30 eigenfields {w\^ 
and 1021,22 are absent), and in this limit all the eigenfields 
coincide or are equivalent to the ones obtained in Ref. Q , 
which are precisely the eigenfields found in the Bona- 
Masso formulation. As emphasized in Ref. *(|, in GR 
taking q = as opposed to s = —1 (which is the value 
proposed in the Bona-Masso formulation), is equivalent 
to the full first order AD MY system. That is, q = 
amounts to not adding the momentum constraints to the 
evolution equations. The ADMY system is then weakly 
hyperbolic. This feature is manifest from Eq. (13.22[) in 
that one has two less independent eigenfields with q = 0, 
so that one cannot reconstruct the original 30 variables 
(a,, di jk ,Kij, V { ) of GR. We then conclude that for STT 
q = also leads to a weakly hyperbolic system since then 
one has only 32 eigenfunctions from which it is impossible 
to reconstruct the 34 field variables listed above. 



gredient to the numerical stability. Nevertheless, it is still 
not very well understood why the conformal approach 
makes such an improvement relative to other strongly 
hyperbolic formulations. More recently the BSSN for- 
mulation has also shown great stability in binary black 
hole simulations when the moving puncture method is 
implemented (24l[25|. Notably, this formulation is able to 
evolve the initial data of two inspiraling black holes until 
the merger stage and beyond without the code crashing. 
This was a long standing problem in numerical relativity. 
Moreover, this formulation has also been very successful 
when matter is included (see for instance (26l. I27I l28l|). 

In view of these advantages, we propose here an analog 
of the BSSN formulation but extended to the STT. More- 
over, we will also deal with the proposed modified Bona- 
Masso slicing condition in this formulation and study the 
roll played by the function in the construction of well 
behaved eigenfields. Just like in Sec. IIIA, we shall de- 
fine new variables so as to obtain a full first order system 
(both in time and space), but following now the confor- 
mal approach of BSSN. 

We consider again the first of Eqs. (|3.4p . and introduce 
the following variables: 





— 6 ^ 7ij ' 


(3.29) 




— P -Ci> a.. 


(3.30) 






(3.31) 


*< 




(3.32) 


3p2 




(3.33) 



B. BSSN Formulation of STT 

The BSSN formulation of GR d 03 is perhaps the 
most popular and useful formulation for numerical ap- 
proaches in use today. This formulation was found to be 
much more stable in numerical simulations than the stan- 
dard ADMY approach. A further analysis of the BSSN 
formulation revealed that this is in fact a strongly hyper- 
bolic system [19(. Although it is clear that the strong 
hyperbolicity property is a key ingredient that causes 
BSSN to perform better than the ADMY formulation 
(which is only weakly hyperbolic), this feature seems to 
be not the only issue in the assessment. In fact there 
are other strongly hyperbolic formulations that perform 
poorly in comparison with BSSN. Apparently it is the 
conformal approach that provides another important in- 



where ip is chosen such that the determinant 7 of the 
conformal metric is equal to one [3!|, 3 f* fe are the 3- 
Christoffel symbols associated with the conformal metric, 
£ is a constant, and Aij is the trace- free part of Kij given 
by: 



Aij = Ktj - i jij K . (3.34) 



From the evolution system of Sec. II one obtains the fol- 
lowing conformal evolution system up to principal part 
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(see the Appendix A for the full system of equations): 



d Qi 
d U 



ad t (sin) 

adiQ 1 , 
-adi 



d a t i 
dadijk 

6 

d K ~ 
doAij ~ - 



f B M K-e^rU 



ae«-WdiA jk 
1 



(3.35) 
(3.36) 

,(3.37) 

(3.38) 
(3.39) 

(3.40) 
(3.41) 



<9 3 r 4 



adi 



(f - 2) e «- 4 )^ i 4i 



Again, we stress that in Eq. (|3.42|) above we have added 
a multiple at; of the momentum constraints (|2.16[) to the 
evolution equation for the variable 3 T\ Notice also A 1 ^ = 
^ %l ^i m Ai m — e( 8 ~^A i: >, and we have also defined 



A 



■Qj) 



f " ' 11 

(3.43) 

where TF means that the quantity between the brackets 
is trace-free (i.e. 7 <J 'A'y = 0). Here c?y := j lm d m ij and 

3 Tj := jji 3 P. For completeness we provide the evolution 
equation for A^-: 



2*, 



f m A l3 + (£-2) S l (i A 



(i " j) 



TF 



BM 



s[id,<)K-$ (0 + 1-0 fya» n i 



(3.44) 

In the above system of evolution equations we consider 
only the reduced system since the lapse and the 3-metric 
propagate along the normal lines. 

Like in the BM hyperbolic formulation of Sec.IIIA, 
there are again 34 variables, namely 15 linearly indepen- 
dent dij m , five linearly independent A\j |4(j, three 3 P, 
one K, three ^ (for the derivatives of the conformal fac- 
tor) , three ai , and finally, four variables II and Qi related 
to the derivatives of the scalar field. 

In order to find the eigenfields associated with these 
34 variables, we consider again a particular direction of 
propagation, say i — x. For q ^ x, it is clear from 
the above evolution equations that there are again 18 



eigenf unctions which propagate with speed — [3 X , namely 
wi,2 = Q q , w 3A = a q , w 5:6 = and u> 7 _i 6 = d qjk . 

Furthermore, we appreciate that the following combi- 
nations 



f 

wn = a x - 6 f B M *x - © y Qx , 



(3.45) 



Wis-20 



3 r + (e - 2) d„ 



r 



(3.46) 



provide another four eigenfunctions with speed — (3 X , 
where d m m1, := -y y 7 d m ij. Thus, we have obtained 20 
eigenfunctions. 

Moreover, just like in Sec.IIIA, we obtain another two 
eigenfunctions from the following combinations 



w 2 i,22 = Q x ± Vt^h 



(3.42) with eigenvalues 



A± = -/3 X T a yff- 



(3.47) 



(3.48) 



Next, combining the evolution equations for K,Q x ,a x 
and II, we form the following two eigenfunctions: 



W 2 3,24 = a 



1 - IBM J J 



e- 1 



1 - Ibm J f 



±e- 2 ^^/f BM l xx K 
with eigenvalues 

Combining the evolution equations for A xx q (A 1 
7* fc A ! fe ■) and A x q , one appreciates that 



(3.49) 



(3.50) 



A xx „±e (c - 2) *i/-7 a « I 



W25-28 



are four more eigenfunctions with eigenvalues 



uj ± = —j3 x ±ae 



-2ip 



(3.51) 



(3.52) 



Furthermore, combining A xxx (A ljk := ^ lk A l \), A x qp , 
A xx and A qp , we obtain the following eigenfunctions: 



W29-32 



A 2 



2j x 



±e «-W^ (a^+JzLa**) >(3 . 53) 



with eigenvalues 



5± = -(3 X ±ae~ 2il 'y/^ xx 



(3.54) 
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As remarked in [6j, W29-32 are in fact only four indepen- 
dent eigenfunctions since these are symmetric and surface 
traceless. 

Finally the last 2 eigenfunctions are obtained by com- 

A xx K 



bining the evolution equations for A 
and II: 

If f> _ 

A XXX Z,XX \ 0~X I J r\X 

^33,34 = A - -7 I la + -jQ 



XXX r,X /")./' AXX 



2£-l 



(3.55) 

where a x := j xl ai, and Q x := j xl Qi. These propagate 
with the following speeds 



?7± = — (3 X ±ae 



-20 



.2£-l 



ryJ. 



(3.56) 



Thus, we have found the complete set of eigenfunctions 
and eigenvalues. The 34 linearly independent eigenficlds 
wi-34 are equivalent to the 34 linearly independent vari- 
ables (Qi, II, Oi, <% fc , if, Aij, 3 f l ). 

Now, like in the Bona-Masso system of Sec.IIIA, there 
are some eigenfunctions (namely ^23,24) that might blow- 
up when /bai = 1. In order to avoid this problem while, 
maintaining the possible choice Jbm = 1, one has two 
simple possibilities: 9 = 1,/bm- These are the same 
conditions found in Sec.IIIA for well behaved eigenfields. 
For instance, with 9 = 1 the eigenfunctions W23.24 reduce 
to 



U>23,24 



L 

f 



Q x ±e' 2 ^y r f, 



BUT 



K 



(3.57) 



with the same eigenvalues, and the blow-up disappears 
even if fsM = 1- The same happens with the choice 9 = 
Sbm and thus we conclude again that in cither case (O = 
1,/bm)j the choice Jbm = 1 requires 9 = 1 which lead, 
as emphasized before, to the modified harmonic slicing 
condition at = -jd*f (cf. Eq. (12~231 with 9 = 

f BM =l). 

The remaining eigenfunctions are smooth for any 
Ibm > and the eigenvalues are real valued. Note, on 
the other hand, that if Jbm = the eigenfunctions 1^23,24 
degenerate and then the system is no longer complete. 

None of the eigenvalues depend on £, and the values 
C = 2, 4 seem the optimal since in that case many ex- 
ponential factors cancel out (C = 4 is the one used in 
the original BSSN formulation of GR, while £ = —2 has 
been used for constructing initial data [18]). On the other 
hand, from Eqs. (|3.55|) and (|3.56[) we appreciate that the 
eigenfunctions and the eigenvalues are real and nondc- 
generate only for £ > 1/2, and therefore the system is 
strongly hyperbolic if this condition holds. A particu- 
lar case corresponds to the value £ = 2 which seems to 



be the most convenient since in that case the eigenfunc- 
tions W18-20 (sec Eq. ({3.460 ). W25-28 (see Eq. (|3.510 ). 
and 1^33,34 (see Eq. Q3.550 ) are simpler. In fact this is the 
value used in the original BSSN formulation of GR. We 
remark that taking £ = corresponds to the case where 
the momentum constraint is not added to the evolution 
equation for 3 f l (cf. Eq. ({3.420 ). and that choice makes 
the system not even weakly hyperbolic since the eigen- 
values (13.560 become imaginary. This is a feature that 
was already present in pure GR @. Actually, when the 
scalar field <p is absent, the hyperbolicity analysis pre- 
sented in this Section is similar (for C = 4 and £ = 2) 
to the corresponding analysis of the standard BSSN for- 
mulation [6j. Finally, it is important to mention that in 
order to avoid any possible divergence in the eigenficlds 
or in the source terms (the terms not contributing to the 
principal part of the evolution equations) of both formu- 
lations, one needs to consider only STT with f{<p) > 
(i.e. a positive definite NMC function). This in turn 
precludes the possibility of having an infinite or negative 
effective gravitational constant G e g (cf. Eqs. (|2.70 ). 



IV. DISCUSSION 

In this paper we have constructed two novel first or- 
der strongly hyperbolic formulations of the STT in the 
Jordan frame along the lines of the BM and BSSN ap- 
proaches. Such constructions show that both formula- 
tions have a well-posed Cauchy problem. This analysis 
fills the gap of a previous study on the Cauchy problem 
of STT [2| and confirms that the Jordan frame is mathe- 
matically adequate for treating the initial value problem. 

One of the most interesting features of the formulations 
presented here is that a modified Bona-Masso slicing con- 
dition is required for the two new systems to be strongly 
hyperbolic while allowing several slicings (/bm > 0) 
which arc natural generalizations of the slicings used in 
pure GR. In particular 9 = 1, /bm in Eq. (|2.240 are two 
simple choices that lead to well behaved eigenfields. 

In the absence of a scalar field, the equations of 
Sec.IIIA,B reduce (for ? = -1, ( = 4 and £ = 2) to 
the known BM and BSSN formulations of GR (when the 
NMC function is trivial, i.e., F(<f>) = 1, the "gravita- 
tional" and the scalar-field sectors decouple completely 
up to principal part). 

What remains to be investigated is the usefulness and 
robustness of these formulations in actual numerical ex- 
periments, as well as the inclusion of a "live shift" . Ac- 
tually, we plan to analyze the dynamical transition to 
the phenomenon of spontaneous scalarization in boson 
stars arising in STT [29( and the subsequent gravitational 
collapse to a black hole with gravitational wave emis- 
sion of scalar type, using one or both of the hyperbolic 
formulations presented here. Both phenomena (sponta- 
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neous scalarization and scalar gravitational waves) are 
even present in spherical symmetry due to the NMC, 
therefore by assuming such a symmetry one can simplify 
the equations without eliminating the interesting physi- 
cal features. 

An important consequence of the analysis presented 
here is that a slightly more general STT which includes 
a function u>(<p) in the kinetic term of the scalar-field 
sector of the action (|2.ip (i.e. the kinetic term has the 
form w(0)(V0) 2 /2) posses a well-posed Cauchy problem 
as well (except for some choices of w(<j>); see Ref. 30}). 
The fact is that all the terms with ll>(4>) do not con- 
tribute to the principal part of the equations associated 
with the metric sector (i.e. the equivalent of Eq. (|2.2p ). 
while the scalar- field sector (the equivalent of Eq. (|2.8p ) 
preserves the quasilinear diagonal hyperbolic form (see 
Ref. [30| for the detailed equations). Thus, up to prin- 
cipal part such STT are identical to the ones analyzed 
here. The relevance of this generalization is that such 
STT can be mapped to the so-called modified theories 
of gravity which are given by a Lagrangian density f(R) 
(R being the Ricci scalar) [3(|. Some specific choices of 
f(R) lead to gravity theories which have been recently 
analyzed in several contexts. Notably, in the cosmolog- 
ical setting such theories have been proposed as an al- 
ternative to dark energy, since they can produce an ac- 
celerating expansion of the Universe without any exotic 
form of matter [UH^]. However, it must be emphasized 
that some of these theories might violate the solar system 
tests [33l. l34j. and some modifications are required to cir- 
cumvent such drawbacks (see Ref. [35[ for a review) . The 
point we want to underline here is that, from the mathe- 
matical point of view, the viability of such theories relies 
heavily on the well-posedncss of the Cauchy problem. 
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V. APPENDIX 

A. Full system of equations 

The full system of evolution equations for the general- 
ized Bona-Masso formulation of STT is the following, 



T t Qi = A («n) , 

4n = a n K + a Q c D c \n (a) + aD c 



fV' - 2f'V 



(5.1) 



4/'(i + 3/")(Q 2 -n 2 ) + i/'T ma tt 



£a = -a 2 (f BU K-e^rU^j 



dt Ul 



a [ Jbm K-e^jH 



-2a Kj, 



(5.2) 
(5.3) 

(5.4) 
(5.5) 



dodijk = -di(a K jk ) + 7 , {j d 2 k)l p l + d m dip 1 

+2d a(3 d k) (3 l , (5.6) 
1% = -DiDj a + a s Rij + aK K. LJ - 2a K a K l } 



+4^G a [ 7 y (S-E)-2 
d V l = -V 1 8il3 l + (1 + q)Di [a (K a - j a K) 

+ \ {d l d m P m - d m 8 m p) 



(5.7) 



aj (K ij - f j K) + 8itG J 1 

i 1 /./ " 1 ,7, ! ■ (5-8) 

where we remind the notations d/dt := dt — and 

d :=dt-^di. 

In fact this last equation can also be written as 

d Vi = Vi d l (3 l + a(l + Oft (k\ - 5\k) 

+ aaj (K{ - 6{K) + \ (d 2 mi f3 m - lin ml d 2 ml P) 

-aK* \2V k ~(2 + ,) 3 T l kl 



(5.9) 



T iifc + (1 + ?) 6 T kij \ - ca8nGoJ., 



where we note that 3 ^ l j k — J l1 (2d(j^)i — dyk) and 
3 T l kl = j^dkji — 9fc(ln^/7). So for the optimal choice 
? = — 1, Eq. (|5.9p reduces to the usual equation of the 
BMSS formulation (cf. Eq. (8) of Ref. [H). 
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We also write here the following terms which contain 
the contributions of the scalar and matter fields @ : 



(S-E) 7ij -2S %3 



Qi Q 3 -(l + /") + f(DiQj + UK tj ) + S? 



7y< 



2/( 



1 



M!2 
2/ 



f /2 



(Q 2 -n 2 )(^-rj-/V 



-2V ( 1 + f — J + (5 matt - £ matt ) ( 1 + y 



47rG a (5 + £) = a f j {D m Q m + U K) 



(5.10) 



n 2 l 



2/ 

^—(1 + 2/") - 
4/ V ^ 2 



3/" 3/" 



2/(1 + ^) 



^matt 



3/ 



/2 N 



-21/ 1 



3/ 



/2 



2/ 



87rGoJi = j 



3/'V* 



-f'(K mi Q m + D t U 



(5.11) 



(5.12) 



The 3-Ricci tensor has the explicit form: 

3 r> n 3W o 3W i 3-riTn 3W 3W 3i^m 

itij — C/ 1 ij — Oi I jl + 1 ij I ml — 1 ira 1 /j 

= -d m d m ij -d (i (2V j} - d j)m m ) 
— A d m m k d^j) 1 + Ad l m (j dj)j m — du m dj m l 
+ 2d l m i (dij m — d m ij) + di m m (%d(ij) 1 — ^(y)) . 

(5.13) 

This expression can be easily rewritten in the form used 
by the BMSS formulation (cf. Eq. (10) of Ref. HH). 



On the other hand, the full system of evolution equa- 
tions for the generalized BSSN formulation of STT is the 



following: 

f t Qi = A (an) , 

4n = aUK + aQ c L> c ln(a) +aD c Q c 



fV' - 2f'V 



(5.14) 



-±/'(l + 3/") (Q 2 -n 2 ) + i/'T ma « 



4<> = ( / BM #-e^n 



/ 



(aK-diP l ) , 



dt^v ~ 2a ^ Aij , 

= -D i D i a + a ( e 2 ^ -4 ^^ A« + i iT 



(5.15) 
(5.16) 

(5.17) 

(5.18) 
(5.19) 



(5.20) 
(5.21) 



-47rG a (-E + S 1 ) 



11. - p-CV 1 



Di Dj a + a ( 3 R i3 - 8vrG S lj 



(5.22) 

TF 



+a 



(5.23) 

where can be extracted from Eq. (|5.10p . In order to 
have a full first order formulation it is understood that 
the second order spatial derivatives of the lapse and the 
physical 3-metric appearing in the above equations have 
to be written respectively in terms of first derivatives of 
di, d i]k and 

The evolution equation for the 3 r* is: 

d 3 T l = a (£ - 2)e^- 4 ^diA u ~a^(^f l diK+ 8wG I' 



-8 l m ai + ^ 3 f % im 



(C + 2)+4-C 



2 ^ 
3 



(5.24) 



13 



where J 1 := 7 J/ whose explicit expression can be ob- 
tained from Eq. (|5.12p . 

We stress that the above equations reduce to the cor- 
responding equations of the original BSSN formulation 
for C = 4 and f = 2 [Il|. 

We must note that for tensor densities T> of weight w 
(here we omit the rank) the Lie derivative has an extra 
weight term: 



«i=0 



wVd m (3 r > 



(5.25) 



where the weight w is related to the power of the deter- 
minant of the physical metric as r y w ^ 2 , and the first term 
on the right-hand side of Eq. (|5.25p is to be understood 
as the ordinary Lie derivative (as though T> were a true 
tensor). For instance, in the case of 7^ and Ay, we see 
from Eqs. (gUP and (j330j) that these are (0,2) tensor 
densities of weight —2/3 and — C/6, respectively so 

(5.26) 

where X>y stands for 7y and Aj with w = —2/3 and 
w = —C/6, respectively. 

For completeness, we provide the expression of the 3- 
Ricci tensor 3 i?y in terms of the conformal metric (|3.29|) 



- + 3R tj 



(5.27) 



where 



2 3 f fc 



T* im a r kh ,(5.28) 



+7 1 i(i 1 j - ) fern 1 ' • <<• 



3 0* 



(5.29) 



where £> stands for the covariant derivative associated 
with the conformal metric 7 13 , and 3 Ffcy := 7/rf 3 r y. In 
fact by putting tildes in Eq. (|5.13[) . and noticing that 

dkrT = 3 y_ l k} = d k (\n^) = (since 7 := 1), we see 
that Eq. (fQg)) can be written as Eq. (|5TT5|) with the 
terms involving d lm m dropped. Moreover, the constraint 
equations (|2.15[) and (|2.16[) can be easily rewritten in 
terms of the quantities defined in Eqs. (|3.29|) — (|3.34|) . 



where the coefficients a — d are in general functions of 
(t,x), and the S UliU2 (ui, 1*2) are source functions. Using 
the judicious guessing approach we shall construct a new 
system of the form (|3.3p which is manifestly hyperbolic, 
and then we shall confront this method with the more 
systematic method which uses linear algebra. 

First we ask for a linear combination w := u\ + o~U2 
to be the eigenfunction. This means that we write an 
evolution equation for w, where cr(a, b, c, d) is a function 
of the coefficients a — d, such that the r.h.s takes the form 
— X a d x w, where A CT will be the eigenvalue. In this way 
we have 



dt (ui + ou-i) ~ —d x (a + ac) u\ + (b + ad) 



, so/ b + ad 

~ - (a + ac) o x u\ H u 2 

\ a + ac 

(5.32) 

Note that with the symbol ~ we have discarded from 
the analysis all the terms which do not contain deriva- 
tives of Ui and u 2 (including the sources). In this way 
we can assume as that the coefficients a — d are "cons- 
tants". By comparing the coefficients of u% on both 
sides of Eq. (|5.32|) . we see that we need to take a = 

d-a±y/(a-d) 2 +4bc 



2c 



(b + ad) 1 '(a + ac), which implies a± 

Moreover, the eigenfunctions will be smooth provided 
c(t, x) does not vanish in some spacetime point. 

Therefore we conclude that the two eigenfunctions 
are w± = u\ + a±ui with the corresponding eigen- 

, \ a+d±V \a-d) 2 +4bc 

values A± = a + a±c = — ; . I he eigenva- 
lues are real and the eigenfunctions are nondegenerate if 
(a — d) + 46c > and c / 0. If this condition holds and 
the eigenfields are smooth then the system is strongly 
hyperbolic [42| . 

Now, let us find the eigenfunctions of the same system 
using the standard method of linear algebra. We rewrite 
the system in a more convenient matrix form as: 



d t 



u 2 



a b 
c d 



d x 



U2 







(5.33) 



The eigenvalues A of the characteristic matrix turn out 
to be A± = a + d± V(° d ) + 4bc ^ w ] ] j c j 1 are identical to the 
ones found before. On the other hand, the eigenvectors 
are: 



B. Simple example of an hyperbolic system 

Let us consider the system 

d t ui + ad x ui + bd x u 2 — S Ul (ui,u 2 ) , (5.30) 
d t u 2 + cd x ui + dd x u 2 — S U2 (ui,u 2 ) , (5.31) 











T 




= 


-0 



v+ = 



Thus, the matrix of eigenvectors reads 



— er_ a + 
1 -1 



(5.34) 



(5.35) 
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Then we can write the eigenfunctions as: 

( fA = w 1 ( Ul )=^—( Ul + a+U3 ), 

(5.36) 

which are th e same found before modulo the factor 
aj \ a = c/ \J (a — d) 2 + Abe. The fact that we do not 
obtain exactly the same eigenfunctions is a consequence 
of the fact that the eigenvectors are unique only up to a 
rescaling. However, both pairs of eigenfunctions w± and 
fi 2 are equivalent. 

Finally, notice that if a = d, then A± = a ± Vbc, 
and w± = Ui ± yfbjc.u%. Thus a plays the same roll as 
— j3 l in the evolution equations of Sec.II. For instance, 
if the coefficients b, c are null, the equations (|5.30[) and 
(|5.31[) decouple and then u\ and ui are themselves eigen- 
functions, indicating that they propagate along the "nor- 



mal lines" (with speed a). This is precisely what hap- 
pens with eigenfunctions Wi-is of Sec.IIIA. On the other 
hand, when b and c are non-null, these coefficients are re- 
lated to the propagation of the eigenfields some of which 
propagate along the "light cones". For instance, if b is 
identified with arf xx (or — qarf 31 , or even afBMj xx ) and 
c with a, in many of the evolution equations of Sec.IIIA 
(assuming propagation only in the l x' direction), then 
the A± are in turn to be identified with the speeds of 
propagation —[3 X ± cty/"f xx or —f3 x ± ay/ '—<;j xx or even 
-f3 x ±ay/f B Ml xx (cf. Eqs. (l3~2Tp . (l3~23l) and (p06]) . res- 
pectively) of several eigenfields (e.g. w± = Ui ± y/"/ xx v.2, 
or w± = ui±y/—qj xx U2, or even w± = u% ± \/ Jbm r ) xx U2 , 
where u\ and i(2 represent field variables like those en- 
tering in the expressions for W23-28, ^29-32 and ^33,34 
respectively) . 
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